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Chiral Ward identities and pion propagation at finite 
temperature in the linear sigma model 



o 
o 
o 

(N 

as: 
< 



> 

00 

o 

O 

o 

a: 

Oh! 



X 



Alejandro Ayala 

Instituto de Ciencias Nucleares, Universidad Nacional Autonoma de Mexico, 
Apartado Postal 70-543, Mexico Distrito Federal 04510, Mexico. 
E-mail: ayala@nuclecu.unam.mx 



Abstract: Working within the linear sigma model at finite temperature, we construct effective 
one-loop vertices and propagators satisfying the chiral Ward identities, for small pion momentum 
and mass compared to the sigma mass. We obtain an expansion for the pion dispersion relation up to 
second order in the parameter m 2 / 'Air 2 f% . At leading order, this expansion reproduces the dispersion 
curve obtained in chiral perturbation theory. We also study pion damping and compute the mean free 
path for pions traveling in the medium before forming a sigma resonance. These results are also in 
good agreement with chiral perturbation theory when using a value for the sigma mass of 600 MeV. 
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1. Introduction 

In recent years, a considerable effort has been de- 
voted to the production and detection of the so 
called quark-gluon plasma -the state of matter 
where the basic QCD degrees of freedom are not 
confined to a single hadron- in relativistic heavy- 
ion collisions. The disentanglement of the signals 
to discern between different outcomes of the colli- 
sion, requires the understanding of the evolution 
of the hadronic state from production to freeze 
out. Since pions are the most copiously produced 
particles in this kind of collisions, a proper under- 
standing of their propagation properties within a 
dense a hot environment represents a key ingredi- 
ent for the understanding of the system's physical 
properties as a whole. 

In general, the hadronic degrees of freedom 
are accounted for by means of effective chiral the- 
ories whose basic ingredient is the fact that pions 
are the Goldston bosons coming from the spon- 
taneous breakdown of chiral symmetry. Chiral 
perturbation theory (ChPT) is one of such theo- 
ries that has been employed to show a number of 
properties exhibited by pions at finite tempera- 



*In collaboration with S. Sahu and M. Napsuciale 



ture, at leading jjj and next to leading ||, [| or- 
ders. A striking result obtained from the second 
order calculations is that the shift in the temper- 
ature dependence of the pion mass is opposite in 
sign and about three times larger in magnitude 
than the first order shift, already at temperatures 
of order of 150 MeV. This result might signal ei- 
ther the breakdown of the perturbative scheme 
at such temperatures or the need to compute be- 
yond next to leading order in ChPT, given the 
large relative corrections found. 

Nevertheless, the simplest realization of chi- 
ral symmetry is still provided by the much stud- 
ied linear sigma model Q . This is a renormaliz- 
able field theory at zero g and (consequently) at 
finite temperature ||. The model has been the 
subject of a renewed interest in view of recent 
theoretical results and analyses of data § 
that seem to confirm the existence of a broad 
scalar resonance, with a mass in the vicinity of 
600 MeV, that can be identified with the cr-meson. 
If confirmed, this last feature could place the 
linear sigma model above the status of a toy 
model to study QCD below the deconfinement 
scale. It thus becomes important to check that 
the results obtained by means of other theories 
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such as ChPT can be reproduced by the linear 
sigma model when working in a given kinemati- 
cal regime. 

Here we show that the linear sigma model 
reproduces the leading order modification to the 
pion mass in a thermal pion medium obtained 
in ChPT at leading order, when use is made of 
a systematic expansion in the ratio m\ / m„ at 
zeroth order -with m (Jl the vacuum sigma 
and pion masses, respectively- and when treating 
the momentum and temperature as small quan- 
tities of order m^. In this regime, the effective 
expansion parameter becomes m 2 /47r 2 / 2 , with 
/tt the pion decay constant. Such parameter al- 
lows for a controlled loop expansion of the pion 
self-energy from where we can study both, the 
modification to the pion mass at next to lead- 
ing order and pion damping in the clastic and 
inelastic channels. This is accomplished by con- 
structing one-loop effective vertices and propa- 
gators satisfying the chiral Ward identities and 
thus respecting chiral symmetry Details of the 
calculation as well as a more extended discussion 
on its implications can be found in Refs. 0. 

2. Effective vertices and chiral Ward 
identities 

The Lagrangian for the linear sigma model, in- 
cluding only the meson degrees of freedom and 
after the explicit inclusion of the chiral symmetry 
breaking term, can be written as j^] 

- AV^+ttV^+tt 2 ) 2 , (2.1) 

where tt and a are the pion and sigma fields, 
respectively, and the coupling A 2 is given by 
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(2.2) 



From Eq. (2.1), we can obtain the Green's func- 



tions and Feynman rules to use in perturbative 
calculations. In particular, the bare pion and 
sigma propagators A„.(P), A a (Q) and the bare 
one-sigma two-pion and four-pion vertices T^ 2 , 



*A 7r (P)<r- 

«A CT (Q) = 



P 2 



-S ij 



lL 12 



iT 



ijkl 



~2i\ 2 {S lJ 5 kl 



S ik 5 jl + 5 il S jk ). (2.3) 



- 04 



are given by (hereafter, capital Roman let- 



These are sufficient to obtain the modification 
to the pion propagator, both at zero and finite 
temperature, at any given perturbative order. 

Alternatively, we can exploit the relations 
that chiral symmetry imposes among different n- 
point Green's functions. These relations, known 
as chiral Ward identities (ChWIs), are a direct 
consequence of the fact that the divergence of 
the axial current may be used as an interpo- 
lating field for the pion ||. Thus, one could 
construct the modification to one of the above 
Green's functions at a given perturbative order 
and from there, build up the induced modifica- 
tion to other Green's functions related to the for- 
mer by a ChWI. For example, two of the ChWIs 
satisfied -order by order in perturbation theory- 
by the functions A 7r (P), A a (Q), T\ J 2 and Y^ 1 
are 

f^Tlf (;0,Pi,P 2 ,P 3 ) = r^(P i; P 2 ,P 3 )^" 

+ if^Ps.Pi)^ 

+ T{ k (P 3 ;P 1 ,P 2 )8 u 
T%(Q;0,P) = f- 1 [A-\Q) A^P)]^, 

(2.4) 

where momentum conservation at the vertices is 
implied, that is P x + P 2 + P 3 = and Q + P = 0. 

At one loop and after renormalization, we 
recall that the sigma propagator is modified by 
finite terms. At zero temperature this modifica- 
tion is purely imaginary and its physical origin 
is that a sigma particle, with a mass larger than 
twice the mass of the pion, is unstable and has 
a (large) non-vanishing width coming from its 
decay channel into two pions. At finite temper- 
ature the modification results in real and imag- 
inary parts. The real part modifies the sigma 
dispersion curve whereas the imaginary part rep- 
resents a temperature dependent contribution to 
the sigma width. 
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In an expansion in the parameter m^/m 2 , 
the leading contribution to the sigma self-energy 
comes from the sigma polarization bubble where 
the intermediate particles are pions. Working in 
the imaginary-time formalism of Thermal Field 
Theory, this is given as 6X 4 f 2 I(uj, q) where the 
function / is defined by 



(uj n - u) 2 + (k - q) 2 + ml 



(2.5) 



Here to — 2mnT and lu u — 2nnT (m, n integers) 
are discrete boson frequencies and q — |q|. From 



Eq. (2Jj) wc obtain the time-ordered version I* 
of the function /, after analytical continuation 
to Minkowski space. The imaginary part of /* is 
given by 
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where Q 2 = q 2 — q 2 , e and are the sign and 
step functions, respectively, and the functions a 
and uj± are given by 



a(Q 2 ) = Jl- 



Am 2 



v Q 2 

q \±a(Q 2 )q 



(2.7) 



whereas the real part of P at Q = is given by 



1 



dk 



Re/ (°) = -^2/ — [1 + 2/(S fe )] , (2. 



where E k = \/k 2 + m 2 , and the function / is the 
Bose-Einstcin distribution 



f(E k ) = 



1 



,E k /T _ 1 
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Therefore, the one-loop effective sigma propaga- 
tor becomes 



*A*(Q) 



ml + 6X 4 f 2 P(Q) 



(2.10) 



The temperature-independent infinities are ab- 
sorbed into the redefinition of the physical masses 
and coupling constants by the introduction of 
suitable counterterms in the usual manner. 

In order to preserve the ChWIs expressed in 



Eq. (2.4), the corresponding one-loop effective 
one-sigma two-pion and four-pion vertices are 

iT\V{Q;P x ,P 2 ) = -2iX 2 f n 6 ij 

[1-3A 2 /*(Q)] , 
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{[1- 3 X 2 I t (P 1 +P 2 )]S ij S kl 
+ [1- 3 X 2 r t (P 1 +P 3 )]S ik 6 jl 
+ [1- 3 X 2 I t (P 1 + P 4 )]5 il S jk }. 

(2.11) 



The functions in Eq. ( 2.11 ) arise from considering 
all of the possible one-loop contributions to the 
one-sigma two-pion and four-pion vertices, when 
maintaining only the zeroth order terms in a sys- 
tematic expansion in the parameter m 2 /m 2 . 

3. Dispersion relation 

We now use the above effective vertices and prop- 
agator to construct the one-loop modification to 
the pion self-energy. Keeping only the leading 
order contribution when considering m^, T and 
P as small compared to m a and to zeroth order 
in ml /ml for which 
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the pion self-energy can be written as 

1 } \2f 2 J 4-7 (2^K 2 +m 2 n 
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P 2 + K 2 



91 



(0)+6I t (P + K)] }, 



(3.2) 



with K = (w„,k) and P — [u>, p). The pion 
dispersion relation is thus obtained from the so- 
lution to 



P 2 + ml + RcIT(P) = , 



(3.3) 
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after the analytical continuation iu> — > po + i ( - 
The temperature dependent infinities contained 
in Eq. ([3^) are exactly canceled by the contri- 
bution from the terms necessary to introduce at 
one loop to carry the (vacuum) renormalization 
procedure |J. 

Let us first look at the dispersion relation at 
leading order, this results from 



P 2 =P 2 + ml [1 + £ g(T/m w )] 



(3.4) 



where £ = m^/Air 2 / 2 <C 1 and the function g is 
given by 

dxx 2 



g(T/m w ) 



\/l + x 2 



1 



This coincides with the result obtained in Ref. JjJ 
by means of ChPT. 

We now look at the next to leading order 



terms in Eq. (3.2). The first of these is purely 
real and represents a constant, second order shift 
to the pion mass squared. The remaining second 
order term shows a non-trivial dependence on P. 
It involves the function S defined by 

d 3 k I\P + K) 



S{P) = T 



(2tt) 3 K 2 + ml 



(3.6) 



The integration in Eq. (3^3) can only be per- 
formed numerically. Including all the terms, the 
dispersion relation up to next to leading order, 
for T ~ and in the small momentum region 
is obtained as the solution to 



2 2 

Po = P 



{l + Zg(T/m x ) + ^-g(T/m n ) 



[9h(T/m v ) - 4g(T/m n )]jml 
e~S{p ,P). (3.7) 



where the function h is given by 
h{T/m^) i 



°° rlh 



(3.8) 



/o E k - 

and S(p ,p) = -(247r 4 )Re5 r (p ,p). Figure 1 
shows the temperature dependence of m T ob- 



tained as the solution to Eq. (3.7) in the limit 



when p goes to zero. Notice that the second or- 
der correction in the parameter £ has the same 
sign as the first order correction. This result is 
opposite to the temperature behavior of the pion 
mass found in Refs. 10, 01. 
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Figure 1: Temperature dependence of m n to sec- 
ond order in the parameter £. For comparison, the 
vacuum value is shown as a horizontal line. 



(3.5) 4. Damping 



Damping is due to two possible kinds of pro- 
cesses: formation of resonances and pion elastic 
scattering. At temperatures on the order of the 
pion mass, pion elastic scattering is subdominant 
since scattering among pions nearly disappears as 
the relative momentum approaches zero. In the 
linear sigma model, this feature is explicitly re- 
alized, for the kincmatical regime discussed here, 
by the description of pion elastic scattering as a 
second order process in the small parameter £. 
We start by first looking at resonance formation. 
This process is the inverse of the one describ- 
ing sigma decay (except that the former does not 
happen at T=0). The relevant diagram to com- 
pute resonance formation is the sunset pion self 
energy with an intermediate sigma line and its 
explicit expression is 



Hi(P) =4.\ 4 f 2 I(P;m„,m 



(4.1) 



where / is the function defined in Eq. (2.5) that 
we rewrite here in terms of the quantities E n = 
\Jk 2 + m 2 and E a = y/(k — p) 2 + m 2 as 

/(P;m - m * )=T ?/||^ 

(4.2) 



K - u>f + El 



The decay rate is given in terms of Imil! by 



rf(?) = 



oPo/T 



ImIIi(po,p) 



( e Po/T _ !) 
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(4.3) 
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where po = \Jp 2 + m\ . Imlli contains both cre- 
ation and decay rates and the factor 



s Po/T /(ePo /T 



1) 



in Eq. (4.3) eliminates the piece describing pion 
creation. We can now compute the mean free 
path A for a pion traveling in the medium before 
forming a sigma resonance. This is given in terms 

by 



of r> 



a = 



(e 



Po/T 



1) 



P 



e Po/T I m rii(po,p) 



,(4-4) 



where v = p/po is the magnitude of the pion's 
velocity. Figure 2 shows plots of A for three dif- 
ferent temperatures and a value of m a = 600 
MeV as a function of the pion momentum. A 
reaches a maximum for p ~ OAm^. Notice that 
the position of the maximum is approximately 
independent of T. Except for the location of the 
maximum, the behavior of the curves is both in 
qualitative and quantitative agreement with the 
corresponding result obtained in ChPT [pi . 




5. Conclusions 

In conclusion, working in the linear sigma model 
at finite temperature, we have found effective 
one-loop sigma propagator and one-sigma two- 
pion and four-pion vertices satisfying the ChWIs 
when maintaining only the zeroth order terms in 
an expansion in the parameter m^/rn 2 ,. We have 
used these objects to compute the two-loop or- 
der correction to the pion propagator in a pion 
medium for small momentum and for T ~ m v 
and showed that the linear sigma model yields 
the same result as ChPT at leading order in the 
parameter £ = m 2 /4ir 2 f%. We have shown that, 
contrary to the result in Refs. & |), the two-loop 
order correction to the pion mass is proportional 
to £ 2 and that this correction is of the same sign 
as the one-loop correction. The shape of the dis- 
persion curve is not significantly altered in the 
kincmatical regime considered. The results for 
pion damping are also in quantitative agreement 
with ChPT when using a value for m a = 600 
MeV. 

Acknowledgments 

Partial support for this work has been received by 
CONACyT Mexico under grant number 32279-E. 



References 



[1] J. Gasser and H. Leutwyler, Phys. Lett. B 184 



(1987) 83; J. L. Goity and H. Leutwyler, Phys 



Lett. B 228 (1989) 517 



p/rr 



Figure 2: Mean free path as a function of momen- 
tum for a pion to travel the medium before forming 
a sigma resonance for m„ — 600 MeV. 



[2] A. Schenk, Phys. Rev. D 47 (1993) 5138 



[3] D. Toublan, Phys. Rev. D 56 (1997) 5629 



[4] J. Schwinger, Ann. Phys. (NY) 2 (1958) 407; M 



Gell-Mann and M. Levy, Nuovo Cim. 16 (1960) 



705 



The pion elastic scattering rate is given in 
terms of the imaginary part of the two-loop pion 



self-energy IT given in Eq. (3.2) 



[5] B. W. Lee, "Chiral Dynamics" (Gordon and 
Breach, 1972). 

[6] L. R. R. Mohan 



Phys. Rev. D 14 (1976) 2670 



lmll(po,p) = -24Vnm6"bo,p) 



[7] N. A. Tornqvist and M . Roos, \Phys. Rev. 

N 



r 2 (po,p) = Imn(p ,p) (4.5) 

Po 

and as anticipated, this rate turns out to be about 



Lett. 76 (1996) 
friep-ph/000813E 



1575 



A. Tornqvist 



[8] M. Svec, \Phys. Rev. D 55 (1997) 5727 . 

[9] A. Ayala, S. Sahu and M. Napsuciale, Phys 



two ordoro of magnitude omallor than Y\_, given 



Lett. B 479 (2000) 156; A. Ayala and S. Sahu 



that it is proportional to £ 2 



Phys. Rev. D 62 (2000) 056007 



5 



